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We performed a search for a light pseudoscalar particle X in the decay if£ — > X, X — > 77 
with the E391a detector at KEK. Such a particle with a mass of 214.3 MeV/c 2 was suggested by the 
HyperCP experiment. We found no evidence for X and set an upper limit on the product branching 



ratio for 



7T°7T°X, X 



77 of 2.4 x 10 at the 90% confidence level. Upper limits on the 



branching ratios in the mass region of X from 194.3 to 219.3 MeV/c J are also presented. 
PACS numbers: 13.25.Es, 12.60.Jv, 14.80.Cp 
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We report the results of a search for the decay — > 
7r°7r°X, X — > 77, where X is a pseudoscalar particle of 
mass in the region 194.3 — 219.3 MeV/c 2 . No experimen- 
tal study for X — > 77 has been published in this mass 
range. This study was motivated by the three events of 
the decay E + — > with the dimuon invariant mass 

around 214.3 MeV/c 2 reported by the HyperCP collabo- 
ration in 2005 Q. 

A sgoldstino interpretation @, d for the HyperCP ob- 
servation showed an upper estimate of the branching ra- 
tio, but this calculation strongly depends on a complex 
coupling constant and only the absolute value can be 
extracted from HyperCP data. The model also showed 
that the X particle's branching ratio (BR) should be 
saturated by the channels X — > 77 and X — > 
with BR(X 77) /BR(X - 10 4 . Another 

model 0] suggested a light pseudoscalar Higgs boson 
interpretation of X, but there was no prediction for 
K° L 7r°7r°X, X 77 decay BR. 

We searched for the K\ — > 7r°7r°X, X — > 77 decay at 
the KEK E391a experiment [5[. Neutral kaons were pro- 



duced by 12 GeV protons incident on a 0.8-cm-diameter 
and 6-cm-long platinum target. The proton intensity was 
typically 2 x 10 12 per spill coming every 4 sec. The neu- 
tral beam [6], with a solid angle of 12.6 /xstr, was defined 
by a series of six sets of collimators and a pair of sweeping 
magnets aligned at a production angle of 4 degrees. A 
7-cm-thick lead block and a 30-cm-thick beryllium block 
were placed between the first and second collimators to 
reduce beam photons and neutrons. The beam size at 
11.8 m downstream of the target, measured with the 
E391a detector, was 3.7 cm (FWHM) including the ef- 
fects of detector resolution. The beam line was kept in 
vacuum at 1 Pa after 5 m downstream of the target and 
1 x 10 -5 Pa inside the fiducial decay region. The K\ 
momentum peaked around 2 GeV/c at the entrance of 
the detector, 11.8 m downstream of the target. 

Figure [T] shows a cross-sectional view of the E391a de- 
tector and defines the origin of our coordinate system. 
The detector components were cylindrically assembled 
along the beam axis. Most were installed inside the vac- 
uum tank to minimize interactions of the particles be- 
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fore detection. The electromagnetic calorimeter, labeled 
"Csl" , measured the energy and position of the photons 
from 7T° and X decays. It consisted of 496 blocks of 
7x7x30 cm 3 undoped Csl crystal and 80 specially shaped 
Csl blocks used in the peripheral region, covering a cir- 
cular area with a 95 cm radius. To allow beam particles 
to pass through, the calorimeter had a 12 x 12 cm 2 beam 
hole at the center. The main barrel (MB) and front bar- 
rel (FB) counters consisted of alternating layers of lead 
and scintillator sheets with total thicknesses of 13.5 X$ 
and 17.5 Xo, respectively, and surrounded the neutral 
beam. Both upstream and downstream ends of MB and 
the upsteam end of FB had Hamamatsu R329-EGP pho- 
tomultiplier tubes [7[ for scintillation light detection. To 
identify charged particles entering the calorimeter, a scin- 
tillation counter (CV) hermetically covered the front of 
the calorimeter. It consisted of a plastic scintillator ho- 
doscope placed 50 cm upstream of the calorimeter with a 
12x12 cm 2 beam hole at the center, and four 6 mm thick 
scintillator plates that connected the beam holes in the 
hodoscope and the calorimeter. Multiple collar-shaped 
photon counters (CC00, CC02-07) were placed along the 
beam axis to detect particles escaping in the beam di- 
rection. The CC02 was a shashlik type lead-scintillator 
sandwich counter with optical fibers running perpendic- 
ularly to the lead and scintillator plates through aligned 
holes, and was located at the upstream end of the K\ 
decay region. The CC03 filled the volume between the 
beam hole and the innermost layers of the Csl blocks 
in the calorimeter. The vacuum region was separated 
by a thin multi-layer film ("membrane") between the 
beam and detector regions. This kept the decay region at 
1 x 10~ 5 Pa despite some outgassing from the detector. 
Further descriptions of the E391a detector are given in 

B|. 

In this analysis, we used data taken in the run period 
from February to April 2005, or Run-II. Data were taken 
with a hardware trigger requiring two or more shower 
clusters in the calorimeter with cluster energy > 60 MeV. 
We also required no activity in the CV and in some other 
photon counters. Because both the decays K\ — > 7r°7r°X, 
X — > 77 and K\ — > 3tt° have a signature of six photons 



in the final state, separation of these decays was crucial 
in this study. 

In the analysis, the K\ -> 3tt° and K\ -> 7r°7r°X de- 
cays were simulated using GEANT3 Monte Carlo (MC) 
simulations Q and were overlaid with accidental events 
taken from the target-monitor accidental trigger. In the 
if£ —> 7r°7r°X decay, X was assumed to decay immedi- 
ately to two photons. To reconstruct K\ — > 7r°7r°X, we 
selected events with six photon-like clusters on the Csl 
calorimeter without any in-time hits on the other detec- 
tors. All the clusters were required to be within the fidu- 
cial region, which was outside the 25 cm by 25 cm square 
around the beamhole and inside a 88 cm circle from the 
center of the beamline. An additional selection criterion 
on the transverse momentum of K\ (Pj- < 0.025 GeV/c) 
was required to suppress missing photon events. After- 
wards, event reconstruction proceeded by solving for the 
decay vertex, assuming the 7r° mass and constraining the 
vertex to lie along the beam axis. The fiducial decay Z- 
vertex (Z vtx ) region was defined to be between 250 and 
550 cm. There were 45 possible combinations to select 
two photon pairs from 6 photons to form two 7r 0, s. The 
most likely pairing was chosen by a minimum which 
was calculated based on the difference between recon- 
structed of the 7r 0, s. In the K\ —> 7r°7r°X event, 
the X mass was the invariant mass of the third photon 
pair (M 56 ) and was reconstructed from the energy and 
hit positions of the remaining two photons and the Z vtx 
determined from two 7r 0, s. The event reconstruction was 
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FIG. 1: (color online) Schematic cross-sectional view of the 
E391a detector. "0m" in the scale corresponds to the entrance 
of the front barrel (FB) detector. K^s entered from the left 
side. 



FIG. 2: (color online) The M56 distribution of the events with 
six photons in the calorimeter. The points with error bars 
show the data, and the red (open) histogram shows the K\ — > 
3tt° MC results normalized to the number of data events. The 
shaded histogram represents the K\ —> 7r°7r°X MC results 
(rax = 214.3 MeV/c 2 ) and was normalized assuming that 
BR (K° L -> 7r°7r°X, X -> 77) is 1.2 x 10~ 4 . 

further optimized by a constrained fit [10] requiring: (a) 
six photons to have the K\ mass, (b) two photon pairs 
to have the 7r° mass, and (c) the if£ momentum vector 
from the target to point to the center of energy of the 
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photons on the calorimeter. The y 2 °f the constrained 
fit (Xfu) was calculated for all 45 combinations, and the 
one with the minimum xju was cnosen to be the correct 
pairing. For successful reconstructions, M 56 equals the 
7T° mass for K\ — >• 37r° decays and the X mass (rax) for 

*2 



7r°7r°X decays. Figure [2] shows M56 of data, the 
ts and the K° L -> 7r°7r°X MC results 
with a 214.3 MeV/c 2 X mass. The signal region is de- 
fined to be 211.3 MeV/c 2 < M 56 < 217.3 MeV/c 2 . Since 
the K\ — > 37T° decay has three 7r 0, s in the final state, 
M56 could be the invariant mass of any one of the three 
7r 0, s. Since only the reconstructed masses of the first and 
second photon pairs are constrained, the pairing with the 
minimum Xfu the worst reconstructed mass of 7r° in 
the third photon pair. This explains the dip in the 7r° 
peak. 

The K\ — > 3tt° mode, with the 19.56% branching 
ratio [Hi, was our dominant background source. Our 



if£ — > 37T° MC results show the tail of the 7r° mass peak 
extended into the X mass region. This background was 
caused by wrong photon pairing combinations, and was 
suppressed by requiring < 4 and x%t < 6. The back- 
ground was further suppressed by rejecting events con- 
sistent with if£ — > 37T° decays. We applied another con- 
strained fit, called the "full constrained fit", by requiring 
the three photon pairs to have ra^o, and reconstructing 
each event as K\ — > 3tt°. As shown in Fig. [3l if£ —> 3tt° 
events in the signal region have smaller x 2 (xjuii) than 
if£ — > 7r°7r°X. By discarding events with Xf u ii < 50, 
the wrong pairing background events in the region M56 
> 165 MeV/c 2 were suppressed by a factor 6.6 x 10 2 . 
The acceptance of K° L 7r°7r°X (rax = 214.3 MeV/c 2 ) 
decays by the Xfuii cu t was estimated to be 61.0% based 
on the simulations. 



With all the selection cuts applied on the data, only 
2 events survived in the 6-MeV/c 2 wide X mass region 
centered on 214.3 MeV/c 2 and 250 events in sideband 
regions (165 MeV/c 2 < M 56 < 211.3 MeV/c 2 and M 56 > 
217.3 MeV/c 2 ). The surviving events had a distribution 
consistent with K\ — > 3tt° decays. We adopted an un- 
binned extended-likelihood method to extract the num- 
ber of signal events. Since the X mass was close to the 
kinematical boundary at 227.7 MeV/c 2 , the background 
shape was modeled by the events in data sideband with 
the ARGUS function [T3 |. which is defined as: 
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exp 



where M ep is the end point of the ARGUS function, x 
is M56, and a determines the curvature of ARGUS func- 
tion. For the signal MC results, the reconstructed mass 
peak is asymmetric because of the kinematical limit of 
the X mass. In this case, a double Gaussian function 
with different mean values of the two composition Gaus- 
sian functions was chosen as the probability density func- 
tion (PDF) to model the signal shape from the signal 
MC results. The modeling result of the 214.3 MeV/c 2 
X mass is shown in Fig. [H The number of signal events 
obtained from the extended-likelihood method, N s , was 
-1.4tJ;g. A dip around 214.3 MeV/c 2 shows that a neg- 
ative number of signal events was yielded. In this paper, 
we searched for X in the mass region > 190 MeV/c 2 and 
no evidence for X was found. We have already published 
an upper limit on the branching ratio of X decaying into 
invisible particles, in the low mass region (rax < 100 
MeV/c 2 ) [ljj. 
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FIG. 3: (color online) The Xfuii distribution of events just 
before the Xfuii cu ^i explained in the text, in the signal region 
of mx = 214.3 MeV/c 2 . The red (open) histogram shows 
the K\ — > 3tt° MC results and shaded histogram shows the 
K° L -> 7T°7T°X MC results (m x = 214.3 MeV/c 2 ). The cut 
point was set at 50. 
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FIG. 4: The modeling result of the 214.3 MeV/c 2 X mass. 
Error bars show the data, and the solid line which is the 
combination of background and signal PDF, is the modeling 
result. The region between two dashed lines is the signal mass 
region. 

Table [J summarizes systematic uncertainties due to (a) 
the signal shape modeling, (b) the signal resolution dif- 



TABLE I: A summary of the systematic uncertainties. The 
uncertainty in the background (bkgd) shape is shown in the 
uncertainty of N s and the others are shown in the percentage 
uncertainty on the branching ratio. 

m x (MeV/c 2 ) 194.3 199.3 204.3 209.3 214.3 219.3 

(a) Sig. shape (%) 6.37 3.39 8.62 9.54 4.07 0.23 

(b) Sig. resolution (%) 1.36 0.29 2.42 2.64 1.79 0.56 

(c) Sig. efficiency (%) 0.93 0.93 0.93 0.93 0.93 0.93 

(d) K° L flux (%) 2.88 2.88 2.88 2.88 2.88 2.88 

(e) Bkgd shape (AN S ) 1.12 1.42 2.35 1.64 0.80 0.27 



The error quoted in S.E.S. was evaluated by summing all 
the systematic uncertainties quadratically. The negative 
central value of the branching ratio was due to the nega- 
tive yield of N s in the signal extraction. The upper limit 
was calculated by integrating up to 90% of the area un- 
der the likelihood function, for N s > 0. The systematic 
uncertainties were incorporated by convolving the likeli- 
hood function with a Gaussian function, and the statis- 
tical uncertainty was incorporated while integrating the 
likelihood function. The upper limit on the branching 



ference between data and the MC results, (c) the signal 
efficiency, (d) the K\ flux estimation, and (e) the back- 
ground shape modeling. The dominant uncertainty was 
source (e). The uncertainty (a) was evaluated by compar- 
ing the value of N s while changing the signal PDF width 
by la fitted error in the signal extraction. The width of 
the 214.3-MeV/c 2 X mass distribution modeled by a sin- 
gle Gaussian function was 1.2 MeV/c 2 . The uncertainties 
(b) and (e) were evaluated in similar way by changing 
the width of the signal PDF and the curvature of the 
background PDF. Since the signal peak was not found 
in data, the uncertainty (b) was evaluated indirectly by 
assuming that the difference in signal resolution is the 
same as the 7r° resolution in reconstructed K\ — > 3tt° 
between data and the MC results. The 7r° mass peak, 
which had a hollow dip in the center of the peak, was 
modeled by the double Gaussian function composed of 
the two Gaussian functions with opposite norm. The 
one with positive norm described the sides of the peak 
and the one with negative norm took out the hollow area 
in the center. The percentage difference of the resolu- 
tion between data and the MC results was estimated to 
be 1.28%. During signal extraction the background PDF 
was fixed. The uncertainty (e) was evaluated by com- 
paring the change in N s while floating the curvature of 
the ARGUS function in the signal extraction. The uncer- 
tainty (c) was simply determined by the statistics of the 
signal MC results. The number of accumulated if£ was 
estimated by the K\ — > 3tt° mode and was cross-checked 
by the K\ — > 7r°7r° mode [11]. Acceptance discrepancies 
in selection cuts between data and MC results create an 
uncertainty in the if£ flux estimation. That uncertainty 
is listed in row (d). The if£ flux at 10 m from the target 
was determined to be (1.32 ± 0.04) x 10 11 based on the 
number of decays downstream of that point. 

The results of different mx are summarized in Ta- 
ble HI N s was the number of signal events yielded, and 
the quoted error was the fitting error. The signal accep- 
tance was calculated using MC code. From the accep- 
tance and the K\ flux, the single event sensitivity for 



7r°7r°X, X — > 77 was defined as: 
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FIG. 5: (color online) The 90% confidence level upper limits 
of the BR(K° L -> 7t°7t°X, X -> 77) for different X lifetimes 
as a function of the X mass. 

ratio for if£ —> 7r°7r°X, X — > 77 decay depends on the 
lifetime (jx) and the mass of X as shown in Fig.0 The 
upper limits were independent of the X lifetime if it was 
shorter than 10 _10 s. The upper limit on the branching 
ratio for K° L 7r°7r°X, X 77 (m x = 214.3 MeV/c 2 ) 
decay was set to be 2.4 x 10 _T in the lifetime region of 
tx < 10 _10 5. The upper limits weaken by a factor of 2-3 
if the X lifetime was 10 _9 s. The upper bound of the X 
lifetime was estimated to be 2.5 x 10 _11 s by [2[. 

We are grateful to the operating crew of the KEK 12- 
GeV proton synchrotron for their successful beam oper- 
ation during the experiment. This work has been partly 
supported by a Grant-in- Aid from the MEXT and JSPS 
in Japan, a grant from National Science Council in Tai- 
wan, from the U.S. Department of Energy and from Ko- 
rea Research Foundation. 



Acceptance • N (if£ decays) 



Present address: Laboratory of Nuclear Problems, Joint 
Institute for Nuclear Research, Dubna, Moscow Region, 
141980 Russia 

Present address: Institute of Particle and Nuclear Stud- 
ies, High Energy Accelerator Research Organization 
(KEK), Tsukuba, Ibaraki, 305-0801 Japan 
Present address: University of Rochester, Rochester, NY 
14627 



5 



TABLE II: A summary of the signal yields (N s ), the acceptances of the signal decay (Accept.), single event sensitivities (S.E.S.), 
central value of the branching ratios (B.R.), and the upper limits at the 90% confidence level(U.L.). 



m x (MeV/c 2 ) 


N s 


Accept. (10" 4 ) 


S.E.S. (10~ 8 ) 


B.R. (10~ 8 ) 


U.L. (10~ 7 ) 


194.3 
199.3 
204.3 
209.3 
214.3 
219.3 


6.4tS;|J 

q c+4.8 
°-°-4.0 

-0 7 +3 4 

u - ' -2.7 

-1 4 +ir 

-0 3 +i r 
u -°-i.o 


1.09 ±0.01 
1.09 ±0.01 
1.09 ±0.01 
1.08 ±0.01 

1.08 ±0.01 

1.09 ±0.01 


7.0 ± 1.3(syst.) 
7.0 ± 2.9(syst.) 
7.0 ± 23.0(syst.) 
7.0 ± 7.6(syst.) 
7.0 ± 4.0(syst.) 
7.0 ± 5.6(syst.) 


44.5+36 4 {stab.) ±8A(syst.) 
24.1+g-j (stat.) ± 10.0(syst.) 
-4.9±li g 8 (stat.) ± 16.3(syst.) 
-10.7t\ 7 2 2 3 (stat.) ± 11.6(syst.) 
-10. Otl^ 8 (stat.) ± b.6(syst.) 

-2A±l* 5 °(stat.) ± 1.9(syst.) 


10.7 
7.9 
5.0 
3.5 
2.4 
2.6 
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